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We study critical point finite-size effects in the case of the susceptibility of a film in which inter- 
actions are characterized by a van der Waals-type power law tail. The geometry is appropriate to 
a slab-like system with two bounding surfaces. Boundary conditions are consistent with surfaces 
that both prefer the same phase in the low temperature, or broken symmetry, state. We take into 
account both interactions within the system and interactions between the constituents of the system 
and the material surrounding it. Specific predictions are made with respect to the behavior of a 
3 He and 4 He films in the vicinity of their respective liquid- vapor critical points. 

PACS numbers: 64.60.-i, 64.60.Fr, 75.40.-s 



I. INTRODUCTION 

Confinement to a finite volume introduces a variety of 
modifications to the behavior of a system in the vicinity 
of a phase transition. The existence of a single bound- 
ing surface leads to surface phase transitions 0, 0, 0, 
critical adsorption Q, , wetting and interface phe- 
nomena . The requirement that the system is of finite 
extent in one or more directions generates effects associ- 
ated with finite-size scaling theory [a, Li [8, 9], including 
shifts in critical points, dimensional crossover, rounding 
of phase transitions and also to such phenomena as cap- 
illary condensation and the interface derealization 
phase transition |lllll2fl . 

In addition, the nature of the interaction within the 
system and between the system and the surrounding 
world influences leading and sub-leading thermodynamic 
behavior at and near a critical point [H Q E M, 13 

In this article we will discuss finite-size effects as they 
apply to the susceptibility of a film of a non-polar fluid 
having a thickness L in which the intrinsic interaction 
J 1 is of the van der Waals type, decaying with distance 
r between the molecules of the fluid as J 1 ~ r -(d+f)_ 
Here d is the dimensionality of the system while a > 2 is 
a parameter characterizing the decay of the interaction. 
The film is bounded by a substrate that interacts with 
the fluid with a similar van der Waals type forces, i.e. of 



the type J 



1,8 



where z is the distance from the 



boundary of the system while a s > 2 characterizes the 
decay of the fluid-substrate potential. For realistic fluids 
d = a = <j s = 3. The discussion in this paper will be 
quite general, but we will be principally interested in an 
Ising type model, which is commonly utilized to represent 
a non-polar fluid. 

As a specific application, we will estimate the param- 
eters required for the prediction of the behavior of the 
finite-size susceptibility in the case of 3 He and 4 He films 
near their respective liquid-vapor critical points. In line 
with an experimental investigation of the phenomena 



that we discuss here, we assume that the film is sur- 
rounded by gold surfaces. We will take into account both 
the van der Waals type interaction between the atoms of 
the 3 He (or 4 He) and the corresponding interaction be- 
tween the Au atoms and 3 He (or 4 He) atoms. 

According to general scaling arguments 0, [13, the 
finite-size behavior of the susceptibility in a film of a 
fluid subject to (+, +) boundary conditions is 




FIG. 1: A schematic {H,T) phase diagram near the critical 
point of a d-dimensional system with (+, +) boundary con- 
ditions [H IU- The line H = 0, T < T c represents the 
bulk phase coexistence line. (In the case of a fluid system 
the vertical -H axis then corresponds to the vapor side of the 
fluid-vapor coexistence curve.) The upper phase coexistence 
line is for a system with thickness Li, while the lower one is 
for a system with thickness L\, where L\ > Li- Away from 
the critical region the shift in the phase boundary is propor- 
tional to I/ -1 , while within it is proportional to L~ A '" . As is 
shown, the lines of a first-order phase transition (with respect 
to the external bulk magnetic field) end at critical points of 
the (d — l)-dimensional system, the positions of which depend 
on L, on the amplitude of the (long-range) surface field and 
on the presence or absence of long-range interactions between 
the spins of the system. 
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where Xbulk(*, ft) is the bulk susceptibility, Xsurface,i(*, ft) 
and Xsurface,2 (t, ft) are the surface susceptibilities, 
£t (T) = Coo(T ™> T+, ft = 0) ~ is the bulk corre- 

lation length, t = (T—T c )/T c is the reduced temperature, 
T c is the bulk critical temperature, £q , a^, ft ftiu,s and 
a w are nonuniversal metric constants, while v, A, r] and 



J 



to are the universal critical exponents for the correspond- 
ing short-range system. For a schematic phase diagram 
in the (H,T) plane of the system see Fig. ffl Assuming 
a > 2 — r\ and a s > (d + 2 — r))/2, one can expand (|l.lfl . 
which leads to 



X(t, ft, L) ~ Xbulk(i, ft) + i 1 [Xsurface,l(*, ft) + Xsurface,2 (* , ft)] 

+Ll /v ^ X ST (a t tL x/v , a h hL A/u - a^L'") + h s L {d+2 - r ' )/2 - a X^attL 1 ^ , a h hL A/l/ ) 
_ L -(a-2 +v ) \hlL-^-^X l l 2 {a t tL l ' v , a h hL A ' v ) + bX^{a t tL l l v , a h hL A ^)\ } , 



(1.2) 



where X s1 ' is the (universal) scaling function character- 
izing the truly short-range system, while the remaining 
part in (|1.2|) describes the contributions due to the (sub- 
leading) long-range part of the interaction. While it is 
well established that X sr tends to zero as exp[— aL/^] 
with a a constant when £/£oo 3* 1, the functions X 1 ^, 
X l * 2 and Xjf are expected to decay only in a power-law- 
m-L way. 

Because of that, whenever L/£oc 3> 1 the functions 
Xf-y , X x * 2 and Xl r will determine the leading-in-L finite- 
size contributions to the susceptibility of the system, and, 
therefore, the leading finite-size behavior of the suscepti- 
bility. Note that due to the lack of "+" "-" symmetry 
in the surface field given (+, +) boundary conditions, one 
finds in the susceptibility a term linearly proportional to 
h w>s (in addition to the term proportional to s ). 

In the case of the "genuine" van der Waals interac- 
tion d = a = a s = 3. For the three-dimensional Ising 
model UK 77 = 0.034, 7 = 1.2385, v = 0.631, a = 0.103, 
(3 = 0.329, 9 = = 0.53. 

In the remainder of this article we will consider only 
the case of identical substrates bounding the system. 
That is, we assume Xsurface,i(i, ft) = Xsurfacc,2(i, ft) = 
Xsnr farafr ft)- Furthermore, for the bulk susceptibility one 
has [23 

Xbulk(£,ft) = Xbulk, sr(i, ft) + Xbulk, lr(*, ft), (1-3) 

where HI El El 



Xbuik,sr(*-+0,ft = 0) =r ± |i|^[l+ (1.4) 

rf \t\° + vf\t\ B ' + vl 2 \t\ 2S + r± a \ty a + ■ ■ ■ ] , 



and 

Xbuik,ir(i - o,ft = 0) = r^i]" 7 [r±|t|^(in|t|)- CT ] . 

_ (1.5) 
Here #2 — 0.90 .25J is the second Wegner exponent, while 
6 a = v{a — 2 + rf) ~ 0.65 is the corresponding exponent 
due to the long-range tail of the van der Waals interac- 
tion. Note that 9 a is only slightly larger than the domi- 
nant Wegner term for the Ising universality class . Apart 
from the logarithmic corrections (of the predicted type 
~ (ln|i|) <T ) the existence of such a correction to scaling 
bulk term has been confirmed within the exactly solvable 
spherical model [13J. The logarithmic term (~ In \t\) has 
been observed, however, only for d + a = 6. In [13 it 
has been demonstrated that such a term is a peculiarity 
of the model which is due to the degeneracy, within this 
model, of the exponents u) = 4 — d and u) a = 6„ jv = a — 2 
(when d + a = 6). Thus, for an Ising like system with 
subleading van der Waals type long-range interactions we 
expect that 



Xbuik,ir(* - o,ft = 0) = r±|t|-T [r±|t| fl *] 



(1.6) 



As far as we are aware, no experimental verification of 
such type bulk correction to scaling has been reported in 
the available literature. 

In the case of (+, +) boundary conditions the surface 
susceptibility is controlled by the extraordinary (or nor- 
mal) surface universality class. One has 



Xsurfaco(i ~> 0, ft = 0) 



r ± 

surface 



t\ 



(1.7) 



with 7s = 7 + 1/ Q. Thus, for the three-dimensional 
Ising universality class j s ~ 1.87. We are not aware of 
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any study of explicit corrections to the behavior of this 
quantity due to van der Waals forces. It is clear that these 
forces do not suffice to change the surface universality 
class. That is, the critical exponents will remain the 
same as in the case of completely short-range interactions 

Let us now investigate in more detail the conditions 
under which the expansion in Eq. (jnj, leading to JOJ, 
can be justified. Some requirements, such as 2— r\— a < 0, 
(d + 2 — T))/2 — a < 0, are obvious and normally are satis- 
fied in any realistic system for which d = a = 3 and r\ <C 1 
(i.e., for the 3d Ising model in which rj ~ 0.034 21] ). Im- 
portant additional conditions arise, however, from the 
fact that we consider a finite system in which power law 
long-range surface fields (i.e. substrate-fluid potentials) 
act. The influence of those long-range surface fields is 
felt everywhere in the finite system, the amplitude of the 
surface field being minimum at the center of the system. 
One can think of this smallest value as a type of a bulk 
field h, which has the effect of displacing the system from 
the position on the phase diagram on which its bulk field 
would otherwise place it. Taking into account contribu- 
tions from both surfaces we obtain for the contribution of 
the long range surface field to an effective bulk symmetry- 
breaking field hh, s = 2h WjS [L/(2£q~)1 . Since the bulk 
magnetic field scales as /i[L/^] A ' i/ one arrives at the 
criterion that in a film the finite-size contributions due 
to the long-range surface fields will be negligible in the 
critical region if 

2\h w , s \ [L/(2£+)p [L/i+] Nv « 1, (1.8) 

i.e. 

2^\h w , s \[m+] A/u -^i. (i.9) 

Note that h w<s > corresponds to attractive walls, 
i.e. walls preferring the liquid phase of the fluid while 
h w .s < corresponds to repulsive walls, i.e. to walls pre- 
ferring the gas phase of the fluid. More detailed discus- 
sion on that point is presented in Appendix |0 where we 
identify h WtS in the framework of a mean-field type model. 
Using the relations between critical exponents it is easy 
to show that Ajv = (d+2 — rj)/2. Thus the relation 
is consistent with the form i|l.l[l . On the other hand, 
A/v = d — (3/v and, therefore, AJv — a = d — o — (3 J v. 
By taking into account that for realistic systems d = a, 
the condition 11. 9|) becomes 

2 ff+1 |/ Vs | [L/C + ]^«l. (1.10) 

For most systems £q" is of the order of 3 A. Taking the 
values of j3 and v to be appropriate to the 3d Ising model, 
we obtain in the case a = 3 

L » Lent = £ + {r +1 \h WtS \Y' P ~ 612 IV.I 1 ' 918 A. 

'1.11: 

Later on in this article we will discuss the determina- 
tion of the magnitude of \h w , s \ for the cases of 3 He and 



4 He films bounded by Au surfaces. For the moment we 
note that in such systems \h w>s \ ~ 4 (see Appendix [UJ. 
The condition L 9000 A must be met if finite-size 
effects due to the van der Waals interactions are to be 
neglible within the critical region of an 3 He or 4 He film. 
If L ~ 15000 A, which is experimentally realizable, then 
one expects van der Waals finite-size effects to play an 
essential role everywhere within the critical region. This 
implies that the value of the finite-size susceptibility at 
T c will depend on L and on the choice of the bounding 
substrate (i.e. on the value of h w>s ). 
Away from T c one expects 

X = Xbulk + "^Xsurfacc + XHamakcr 1 ' ' (1-12) 

In a fluid system below T c the Hamaker term com- 
prises contributions due to both the fluid-fluid and the 
substrate-fluid interaction (the substrate-substrate inter- 
action does not contribute because the substrate density 
p s does not depend on the magnetic field) . In a fluid sys- 
tem the density is the quantity that couples via the van 
der Waals type interaction. The average local density is 
non-zero both below and above the critical point, whereas 
in a magnetic system mo(T > T C ,H = 0) = 0. Because 
of this, there exists a difference in the finite-size behavior 
of the fluid and the magnetic systems when T > T c . One 
can argue that for a magnetic system in which T > T c 
and h = [H 

1 1 , 

X = Xbulk + ^Xsurface + J^+lX H U-13) 

Below T c Eq. (| 1 . 1 2|) is again valid. As we see, the terms of 
the order of L~ ( - d+ ^ will never be important if substrate 
potentials are present in the system. Such terms, how- 
ever, give rise to the leading finite-size contributions away 
from the critical region in the case of periodic boundary 
conditions 0, ITS Il5l ITsj . Contributions from the fluid- 
fluid potential are proportional to m 2 and, therefore, will 
be of importance only below the critical temperature (if 
\H\ <C 1). They then constitute one of the Hamaker 
terms. 

Expressions similar to Eq. (|1.2|) have been shown 
to describe the finite-size behavior of the susceptibil- 
ity in a fully finite system subject to periodic boundary 
conditions — in which case h W)S = — both in the instance 
of the exactly solvable spherical model [l3LllS . ll^| and via 
e-expansion techniques (up to first order in e), in 0(n) 
models [lq . 

Figure |3 displays the behavior of the susceptibility in 
a film in which the interactions are completely short- 
range and of a film in which both the fluid-fluid and the 
fluid-substrate interactions are long-range, the latter case 
corresponding to the actual experimental situation. We 
observe that the curves behave quantitatively differently 
in the two systems. At T c and at coexistence, the sus- 
ceptibility of the van der Waals type system is severely 
suppressed in comparison to that of the short-range sys- 
tem. As explained above and as we will see in more detail 
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below, the magnitude of the maximum depends on the 
strength of the substrate-fluid coupling and is, therefore, 
not universal. Furthermore, there is a shift in the po- 
sition of the maximum of the susceptibility. We expect 
that the curves shown here will resemble those obtained 
via experimental investigations. 
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FIG. 2: A comparison between the behavior of the scaling 
functions of the susceptibility X x (xt) for a fluid in which 
both the substrate-fluid and the fluid-fluid interactions are 
short-range and the same quantity for a fluid in which both 
interactions are long-range. The thickness of the film is La, 
where L = 3000, while a is the average distance between the 
molecules of the fluid. The scaling variable is xt = tL 1 '", 
where r = 1 — T c /T with T c being the critical temperature 
of the corresponding short- or long-range system. The curves 
differ essentially from each other, in that they have different 
magnitudes at T c . The maximum of the scaled susceptibility 
is also shifted. 
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In this article, we will explore the mechanisms under- 
lying the similarity of the susceptibility in a fluid film 
near its bulk critical point to the susceptibilities plotted 
in Fig. |2 We investigate the behavior of the suscepti- 
bility as a function of the temperature, the external bulk 
magnetic field, the thickness of the film and the strength 
of the fluid-substrate interaction. 



The structure of the article is as follows. First, in 
Section [n] we present a precise formulation of the model 
of interest along with the analytical expressions needed 
for its numerical treatment. The principal results for 
the behavior of the finite-size susceptibility are presented 
in Section IIIII Some technical details — including sums 
containing van der Waals type potentials — can be found 
in Appendix^ An estimation of the parameters needed 
in order to apply the current results for the behavior of 
the susceptibility in 3 He and 4 He films is contained in 
Appendix IO 

II. THE MODEL 



Perhaps the simplest model that captures the basic 
features of systems with van der Waals interactions is a 
modification of the one utilized by Fisher and Nakanishi 
in their mean-field investigation of short-range systems 
|27l l28j . One starts with the following form of the func- 
tional for the free energy of the lattice system 



1 + m(r) 



In 



1 + m(r) 



1 — m(r) 



In 



1 — m(r) 



| — h(r)m(r) — — K(r, r')m(r)m(r'), 



(2.1) 



r 



where K{v,r') = /3J (r, r') is the non-local coupling be- 
tween magnetic degrees of freedom, h(r) is an exter- 
nal magnetic field and the magnetization m(r) is to be 
treated as a variational parameter. Note that in Eq. (|2.1|) 
the term in curly brackets corresponds to the entropic 
contributions, while the other terms are directly related 
to the interactions present in the system. 

The variation of (|2.1|) with respect to m(r) leads to the 
following equation of state for our system 



m*(r) = tanh 



^K(r,r')m*(r') + /i(r) 



(2.2) 



This equation can be solved numerically by applying the 
Newton-Kantorovich method. One is able to treat rea- 



sonably thick films, with L/a of the order of 3000 layers, 
corresponding to the experimental setup that we envis- 
age as an example in our study. Its solution for a given 
geometry and external fields h(r) determine the order- 
parameter profile m*(r) in the system. 



We will be interested in a system with a film geom- 
etry. Then if r = {rii,^} and h(v) = h(r»,z) — h(z) 
one has, because of the symmetry of the system, m(r) = 
m(r\\ , z) = m(z). The magnetization profile now depends 
only on the coordinate perpendicular to the plates bound- 
ing the van der Waals system. In this case, Eq. Q2.2[l . 
which describes the behavior of the magnetization profile, 
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becomes 



m*(z) = tanh 



/3^J(r -r')m*{z') + h(z) 



, (2-3) 



where r' = (rj.,2:'). Obviously, the above equation is 
equivalent to 



m*(z) = tanh 
where 



PJ2 J{z-z')m*{z') + h{z) 

. z'=0 



j{z) = y j j{h-^\^) = Y. j ^p 2 )- 



(2.4) 



(2.5) 



We will now assume that the fluid molecules occupy 
the region in space characterized by < z < L and that 
the layers z — and z = L satisfy the (+, +) boundary 
conditions, i.e. m(0) = m(L) = 1. The number of layers 
containing spins that can fluctuate is, therefore, L — 1. 

Equation (|2.4|> is equivalent to 

arctanh[m*(z)] = h{z) + (3^J {z - z')m* {z 1 ). (2.6) 

z' 

Taking the functional derivative of the both sides of Eq. 
(12. 6|) with respect to the field h{z*) applied to the layer 
z* , we obtain 



G(z,z*) 

1 II \ ~~ °z,z* 

1 — m yz) 



/3^i(z-z')G(z',z*), (2.7) 



where 



G{z,z*) 



(2.8) 



5m(z) 
Sh(z*) 

J2{S(0,z)S(rlz*))-{S(0,z)){S(rl,z*)) 



Equation i|2.7[l can be rewritten in the form 

Sz.z' 



E 



1 — m 2 (z') 



(3J{z - z 1 ) 



G{z',z*)=5 z , z *. (2.9) 



Now it is clear that the solution of the above equation 
with respect to G(z', z*) is 



G(z',z*) = (R- 1 )„, 



(2.10) 



where R 1 is the inverse matrix of the matrix R with 
elements 



R z 



i^(z') 



(2.11) 



1 — m 

One can define the "local" susceptibility 

Z* 

= *£(S(0,z)S(rl,z*))-(S(0,z))(S(rlz*)). 

r* 

(2.12) 

Obviously, \ = Tl z x( z )/ (L + l) is the total susceptibility 
of the system per unit spin (see Eq. [j2.12J) ). Thus, one 
has 



X (z)=£(R- 



and 



_|_ \ ^ / )z,z* 



(2.13) 



(2.14) 



In Appendix lAl we demonstrate that the function j{z) 
can be written in the form 

J(z) = J 1 [cd-iS(z) + c™ ! [5(z — 1) + S(z + 1)] 

+g d (z)6(z-2)}, (2.15) 

where S(z) is the discrete delta function, while 9{z) is the 
Heaviside function. Explicitly for d = a = 3 one has (see 
Appendix |A"|) 



C2 



y — 

^ 1 + |n|' 



3.602, 



(2.16) 



nn 
c 2 



-— 7T 

3 



(-l) 1/3 K (^2-2iV3ir) - (-l) 2/3 K Q (\/2 + 2iV3ir) 



and 



r V3 , f 1 \ 1 , / I I 
V 3 arctan — ^ — In I 1 H — - I + - m | 1 — 



(2.18) 



The layer magnetic field h(z) is the only quantity in Eq. We take it to be of the form 
(|2.6|l the exact form of which still requires specification. 
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where h WjS reflects the relative strength of the fluid- wall 
and fluid- fluid interactions (see Eq. (|(J12|) below). The 
above expression is consistent with the fact that the sub- 
strate occupies the region R d ^ 1 x [L + l,oo] U K d_1 x 
[— (L+l), — oo]. For 3 He and 4 He bounded by Au surfaces 
we show in Appendix IH1 that h ws — 4. 



III. NUMERICAL RESULTS AND FINITE-SIZE 
SCALING ANALYSIS 

To determine the total susceptibility \(T, H\L, h WyS ) 
and its "scaling function" 



X x = L-^ v X (T,H\L,h w<s ) 



(3.1) 



in a fluid film with thickness L, one has to solve Eq. 
()2.6(l for 1 < z < L— 1, which allows one to construct the 
matrix R with the use of Eq. (|2.11|) . Finally, on the basis 
of Eq. I|2.13[l one obtains the local susceptibility x(z), 
1 < z < L — 1, and, summing, the total susceptibility. 
We recall that within the mean field treatment one has 
/3 = v = 1/2, 7 = 1, i] = and, for the extraordinary 
surface transition, 7 S = 3/2 |l|. 

The analytic solution to the set of coupled nonlinear 
equations for the magnetization profile is, at present, 
known only for purely short-range, continuous, systems 
[2^ | . We review this solution in section llll Al Even in this 
case the the finite-size susceptibility has to be determined 
numerically. The results are, again, presented in section 
1111 Al Numerical methods appear to be unavoidable in 
order to solve the equations for the magnetization profile 
in the case of the long-range interactions of the van der 
Waals type. The results in this case will be discussed in 
section llll Bl 



A. The model with purely short-range interactions 
for H = 



The equations to be solved in the continuum version of 
the purely short-range model of a mean-field Ising strip 
under (+, +) boundary conditions are 



1 







1 







dz 2 
dz 2 



m 2 (z) 
—m 2 (z) 



X(z), (3.2) 
m[z). (3.3) 



Because conditions are identical at both bounding sur- 
faces of the system, the solutions of the above equation 
satisfy m'(L/2) — and y'(L/2) = 0. The magnetization 
profile is known exactly 29] : 



a) when tL 2 > ir 2 [with t = (T - T c ) /T c ] 



m(z) 



2K(k) dn(C;fc) 
L sn(£; k) 



(3.4) 



where 

tL 2 = [2K(k)] 2 (2k 2 - 1), C= [2K(k)/L]z, (3.5) 

and k 2 > 0. 
b) when tL 2 < it 2 

2K{k) 1 



m(z) — 



L sn(C; k) ' 



(3.6) 



where 



tL 2 = -[2K(k)] 2 (k 2 + 1), 
and k 2 > 0. 



{ = [2K{k)/L]z, (3.7) 



Here K(k) is the complete elliptic integral of the first 
kind, dn(£; k) and sn(£; k) are the Jacobian A amplitude 
and the sine amplitude functions, respectively. The bulk 
critical point T — T c corresponds to k 2 = 1/2. The above 
expressions are consistent with the following scaling form 
for the order parameter: 



m( 2 )=i-^X m ( I ,tL 1 ^), 



(3.8) 



with (3 = v = 1/2. 

The results from the finite-size scaling analysis of 
the behavior of the susceptibility for a system with 
short-range interactions are summarized in Figs. |3| 
and Fig. [3] compares the finite-size susceptibilities 
L~f/" X (T, h\L, h w>s ) for short-range films with L = 100 
and L = 3000 layers. Fig. 0] presents the corresponding 
results for the ratio x(T,h\L,h WyS — 0)/x(T — T c ,h — 
0\L, h WiS = 0). In both plots the scaling variable is 
Xt = (1 — T c /T) L x l v . The curves demonstrate a reason- 
ably good scaling where the small deviations for L = 100 
from the L — 3000 curve can be explained with the am- 
biguousness in the definition of L (for a discussion see 
Appendix [BJ), as well as with the corrections to scal- 
ing terms and the role of the background (nonuniversal 
terms). Note that the curves present the behavior of the 
total susceptibility and not only of its singular part. 



B. The model with van der Waals type interactions 

We note that the critical temperature T c depends on 
the presence or absence of long-range fluid-fluid interac- 
tions in the system. Let us denote by T Cj i r the critical 
temperature of the system with subleading long-range 
fluid-fluid interactions (from van der Waals type) and 
with T c sr the corresponding temperature for short-range 
fluid-fluid interactions. If K — (3 J 1 then it is straightfor- 
ward to show that K c \ r ~ 0.161, while i^ c , sr ~ 0.168 in 
the framework of our model defined by Eq. (|2.1|) . 

The results from the finite-size scaling analysis of the 
behavior of the susceptibility for a system with long- 
range van der Waals type interactions are presented in 
Figs. El El El El l^landlTUI The scaling procedure is ex- 
plained in details in Appendix iBl 
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FIG. 3: The upper dotted line, for L = 100, and the lower 
dotted line, for L = 3000, illustrate the behavior of the scal- 
ing function of the susceptibility X x (x t ) for a fluid in which 
both the substrate-fluid and the fluid-fluid interactions are 
short-range. The scaling variable is x t = (1 — T c /T) L 1 ^ . 
Note that the maximum of the finite-size curves is at T > T c 
which is due to the stabilizing effect of the (+,+) boundary 
conditions on the order parameter in the finite system. We 
observe that the both curves differ at most around their max- 
imal value. As it will be demonstrated in Fig. 2]the principal 
reason for this deviation is the improper choice of the value 
of the distance L between the plates. 



1. The temperature dependence at H — 

Let us first consider the behavior of the finite-size sus- 
ceptibility in the absence of an external magnetic field, 
which is equivalent to the behavior in a fluid system along 
the liquid-vapor coexistence line. Our results are pre- 
sented in Figs. |S]and[f)l One observes that the suscepti- 
bility possesses a maximum above the bulk critical tem- 
perature (which reflects that fact that the (+, +) bound- 
ary condition stabilize the long-range order at tempera- 
tures slightly above T c ) and that this maximum is weaker 
than for the corresponding short-range system (see Fig. 
©). The scaling variable is x t = (1 - T c /T) L 1 ^. Note 
that the scaling functions decay much more slowly as a 
function of \xt\ in comparison with the short-range case. 
The maximum of the short-range case is around xt — 4 
while in the case of a van der Waals type system it is 
around x% ~ 8. These results imply that, as expected, 
the long-range tails of the interactions help to stabilize 
the long-range order even a bit above the corresponding 
limit in T for the short-range system with (+, +) bound- 
ary conditions. 



2. The temperature dependence of the susceptibility at 

In this section we consider the behavior of the finite- 
size susceptibility for values of the external bulk magnetic 



FIG. 4: The lines show the behavior of the scaling function 
of the normalized susceptibility X x (xt)/X x (0) for a fluid in 
which both the substrate-fluid and the fluid-fluid interactions 
are short-range. The curves are for L = 100 and for L — 3000. 
The scaling variable is x t = (l-T /T)L 1,v . Both curves 
coincide perfectly well with each other. The main reason for 
the deviation of the curves in Fig. El is the improper choice of 
the length L (for a discussion see Appendix |5J. 



field that support the vapour phase of the fluid. The 
results are presented in Fig. Note that by changing the 
sign of H (by choosing negative values of Xh = (3HL A ^ V ) 
one can show that the position of the maximum in the 
behavior of the susceptibility, which for H — is at T > 
T c , moves gradually toward T c and for negative field with 
large enough magnitude can even be below T c . In Fig. 
we present the above results as a function of L/^ t - This 
resolves the question about the value of the nonuniversal 
metric factor in x t . Note that then the shape of the curve 
and the position of its maximum shall be a reasonable 
approximation for the real experimental system of 3 He 
or 4 Hc films bounded by Au surfaces. 



3. The field dependence of the susceptibility at T 7^ T c 

Here we analyze the behavior of the finite-size suscep- 
tibility as a function of the field scaling variable Xh = 
(3HL A / U for fixed values of the temperature close to the 
bulk critical temperature T c . The results are presented 
in Fig. M The curves are for r = 0, 1(T 6 , 1CT 5 , 1CT 4 . 
Note that the behavior of the susceptibility possesses a 
peak that essentially differs from the corresponding back- 
ground contribution only for |r| < 10~ 6 . In Fig. 1101 we 
present the same results, this time as a function of L/£h- 
This resolves the question about the value of the possible 
nonuniversal metric factor in Xh- Thus, the shape of the 
curve and the position of its maximum shall be a reason- 
able approximation for the real experimental system of 
3 He or 4 He films bounded by Au surfaces. 
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FIG. 5: The upper dot line (for L = 100) and the lower dot 
line for (L = 3000) show the behavior of the scaling func- 
tion of the susceptibility X x (xt) for a fluid in which both 
the substrate-fluid and the fluid-fluid interactions are of van 
der Waals type, i.e. are long-range. The scaling variable 
is x t = (1 - T c /T) L 1/v . We observe that both curves dif- 
fer at most around their maximal value. The amplitude 
of the surface field for L = 100 is h W:S = 0.73 while for 
L — 3000 it is h w s — 4 which ensures that the variable 
2^+ 1 |^, s |[L/C +] A/ ^ <Ts = 2 4 ' 5 |/ W |[L/£+]- 1/2 has the same 
value for both the cases. 



FIG. 6: The lines show the behavior of the scaling function 
of the normalized susceptibility X x (xt)/X x (0) for a fluid in 
which both the substrate-fluid and the fluid-fluid interactions 
are long-range. The curves are for L — 100 and for L — 3000. 
The scaling variable is xt = (1 — T c /T) L ly/ " . We observe that 
the both curves coincide perfectly well which each other only 
near the bulk critical point. The deviation of the curves from 
each other is due to the effect of the van der Waals interaction 
(compare with the short-range case) and increases with the 
increase of \xt\. Note also that the scaling function decays 
much slower with \x t \ in comparison with the short-range case. 



IV. CONCLUDING REMARKS AND 
DISCUSSION 

In this article we investigated the behavior of the sus- 
ceptibility in thin films with van der Waals type long- 
range interactions. Prominent examples of such systems 
are simple nonpolar fluids in thermodynamic equilibrium 
with their vapor, as well as binary fluid mixtures close to 
their demixing point. We have seen that despite the fact 
that this kind of interaction does not change the critical 
exponents of the system it nevertheless gives rise to a 
variety of finite-size effects that become dominant when 
L/^oc ^ 1- Furthermore, we have formulated a criterion 
regarding the conditions under which such effects also 
essentially modify the finite-size behavior of the suscep- 
tibility everywhere within the critical region. According 
to this criterion if the thickness L of the film is such that 
(see Eq. (ED)) 

L<L crit = z+(m\h w , s \y^, (4.i) 

the effects due to the van der Waals interaction are es- 
sential and cannot be neglected. Here £q" is the ampli- 
tude of the bulk correlation length while the (dimension- 
less) factor h Wt s characterizes the strength of the fluid- 
substrate interaction ~ h WtS z~ a ', where z is the distance 
from the substrate surface toward to the bulk of the 
fluid. In the case of 3 He or He bounded by Au sur- 
faces we find that at the corresponding critical point of 
both the fluids h w<s ~ 4 (see Appendix . Then, for 
a three-dimensional Ising type system the direct evalu- 



ation yields the estimate L C rit — 9000 A. A comparison 
when L/a = 3000 layers (where a is the lattice constant) 
between the behavior of the finite-size susceptibility of 
a system with completely short-range interactions and 
one with long-range fluid-fluid and substrate-fluid inter- 
actions is given in Fig. [21 One observes a clear distinc- 
tion between the curves in the critical region. The cal- 
culations were preformed using a mean-field type model 
which is described in detail in Section ^ In order to 
determine the susceptibility one solves L/a — 3000 cou- 
pled nonlinear equations; we make use of the Newton- 
Kantorovich method. We have chosen a film with this 
number of layers because it corresponds to a realistic ex- 
perimental setup of 3 He film between the Au electrodes 
of an experimental cell in which the smallest distance be- 
tween the plates is 1.5 /im; the distance r$ between the 
atoms of 3 He at its liquid-vapor critical point has been 
estimated to be r ~ 4.9 A (see Appendix EJ. For such 
a system the behavior of the susceptibility as a function 
of (1 - T c /T)L l ' u for different values of the bulk field 
x h = [3HL A / U is shown in Fig. □ while in Fig. the 
susceptibility is plotted as a function of Xh for a few fixed 
values of t. The same data are shown in Figs. I%1 andHTH 
as functions of L/£t an d L/^h, respectively. We expect 
these curves to resemble the actual experimental data 
for 3 He (or 4 He) film. One observes that for Xh = 0, 
the maximum of the susceptibility is above the bulk crit- 
ical point. Furthermore, we note that when one is on 
the vapor side of the bulk phase diagram the maximum 
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FIG. 7: The lines show the behavior of the scaling function 
of the susceptibility X x (xt, Xh) for a fluid in which both the 
substrate-fluid and the fluid-fluid interactions are long-range. 
The substrate-fluid potential is characterized by h w . s = 4, 
which correspond to the situation of 3 He or 4 He films bounded 
by Au surfaces (see Appendix[njl- The curves are for L = 3000 
and at x h = 0,-5,-25,-50,-100, where x h = /3HL A/v . 
The scaling variable is x t = (1 - T c /T) L 1/v . The changed 
shape of the curve for Xh = —100 signals that this curves is 
the precursor of a first-order phase transition. 
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FIG. 8: The same as in Fig. |7| but the scaling function is 
presented as a function of L/£t at several fixed values of 
xh- We expect this curve to be a good approximation of 
the corresponding finite-size scaling function for a real three- 
dimensional experimental system of 3 He or 4 He films bounded 
by Au surfaces. 



moves toward the bulk critical point with increasing dis- 
tance from the coexistence curve in the phase diagram 
as a function of —Xh and that for — Xh large enough the 
maximum positions itself at a temperature less than T c . 
Additionally, one observes that as a function of the scal- 
ing field variable {xh, or = L/^h) the susceptibility 
does not display an easily distinguishable maximum un- 
less \t\ < 10 -6 . The position of this maximum changes 



FIG. 9: The lines show the behavior of the scaling function 
of the susceptibility X x (xt,Xh) as a function of Xh at several 
fixed values of x t for a fluid in which both the substrate-fluid 
and the fluid-fluid interactions are long-range. The substrate- 
fluid potential is characterized by h w , s = 4, which correspond 
to the situation of 3 He or 4 He films bounded by Au surfaces 
(see Appendix ICl . The scaling variable is Xh = flHL A ^ . The 
numerical calculations are performed for L — 3000 layers. 
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FIG. 10: The same as in Fig. but the scaling function 
is presented as a function of L/£h at several fixed values 
of r. We expect this curve to be a good approximation of 
the corresponding finite-size scaling function for a real three- 
dimensional experimental system of 3 He or 4 He films bounded 
by Au surfaces. 



very slightly with t and is around sign(h)L/£h — —8.7 on 
the vapor side of the bulk phase diagram. 

We stress that all of the above curves depend on the 
value of h WjS and are, thus, nouniniversal. This implies 
that if h WjS is kept fixed when L changes one will obtain 
different curves for the different L's. The same will also 
be true when L is kept fixed but h WtS changes. This is 
illustrated in Appendix iBl — see Figs. 1121 and IT51 Only 
when h W)S /y/L is kept fixed can a curve that does not 
depend on L nor on h w s , but just on their ratio, be ob- 
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tained - see Fig. In practice, one does not know the 
precise value of the system size L - an important issue 
for very thin film thicknesses (for more details see Ap- 
pendix [0). When this is also taken into account further 
improvement of the data collapse can be achieved as seen 
in Fig. 1151 where the corresponding data for the suscep- 
tibility for fixed h WjS /*/L are normalized by its value at 
the bulk critical point. Definitely the above predictions 
are clearly experimentally verifiable. For a given fluid, 
one can either change the size L of the distance between 
the plates of the experimental cell (i.e. the fluid film) 
or the corresponding substrate (i.e. the value of h w _ s ). 
The bigger the change in the new h WjS , the bigger will be 
the deviation of the new curve of the finite-size suscepti- 
bility from the old one. Of course, this will be possible 
only if L < £ C rit- When L 3> £ C rit the universal behav- 
ior of the finite-size susceptibility in its standard form 
will prevail and the van der Waals interaction will only 
lead to small, probably experimentally negligible, correc- 
tions to the universal curve. Let us also note that when 
the experiment is performed for 3 He or 4 He around their 
liquid-vapor critical points according to our predictions 
one will obtain for any fixed L - even when L < L cl -i t 
- practically the same curve X x for the finite-size sus- 
ceptibility of both the fluids. For L < L cl -it this will be 
simply due to the fact that h WtS ~ 4 for both the flu- 
ids (see Appendix Q , while for L ^> L cr j t that will be 
due to universality since near their liquid vapor critical 
points these fluids belong to the three-dimensional Ising 
universality class. 

In the work described above, the effects of retardation 
on the van der Waals force have been neglected. These 
effects set in for distances r larger that 160 A between 
the He and Au atoms [3(| • A possible concern is the in- 
fluence of retardation on the system investigated in this 
paper. We have estimated that influence by performing 
numerical calculations in which an "extreme" retarda- 
tion was imposed, in that the interaction potential was 
set to zero for separations r > 100 layers. We discovered 
that the numerical consequence of retardation grows with 
increasing film thickness. At L = 3000 the difference be- 
tween the calculation with the unretarded van der Waals 
force and the force with long-distance cutoff is 13.5% at 
T = T c . This is an overestimate of the actual influence 
of retardation in experimental realizations. We note that 
retardation reduces the suppression of the susceptibility 
by the van der Waals interaction — see Fig. |5| 

Finally, let us note that in the current study we did not 
take into account the influence of gravity on the behav- 
ior of the finite-size susceptibility. Because gravity leads 
to stratification of the density one expects an additional, 
more complicated z dependence of the local and, there- 
fore, of the total susceptibility. We hope to return to this 



question in a future publication. 
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APPENDIX A: DERIVATION OF THE 
EQUATION FOR THE MAGNETIZATION 
PROFILE IN A VAN DER WAALS FILM 

Here we derive an explicit form for the function J{z) 
for some basic cases of special physical interest. Namely 
we take J(r) to be of the "van der Walls form" 



J(r) = 



J 



^ _|_ r d+o ' 



(Al) 



where r = |r|, and d = a = 3 for the "real" (nonretarded) 
van der Waals interaction. 

Then one can further simplify the sum in the right- 
hand side of the above equation. Using the identity 

1 f°° 

-—-= dt cM-zt)^ 1 E a , a (-t a ), (A2) 
l + z J 



where 



E a , f ] (z) = 



k=0 



T(ak + f3)' 



a > 0, (A3) 



are the Mittag-Leffler functions, the sum can be rewritten 
in the form 



£j(r-r')m(z') = J^z - z')m(z'), (A4) 



z'=0 



where 



Jd,a(z) 

= J dt t^)' 2 - 1 E d+ „ d+ ,{-t^/ 2 ) 
Jo 2 ' 2 



E 



e~ tv \ e' tz 



(A5) 



The main advantage of the above form is that it factorizes 
the summation over the components of rji . With the help 
of the Poisson identity Eq. (|A4|) becomes 
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£j(r-r')m(z') = 



m{z) 



(A6) 



r , , /7r\( d -!)/2 d+ „ j 

/o ^ * 



r^" 1 £4+^ (-^) E e-*(*-*') a m(z / ) 



, /7T\(rf-l)/2 

I i 2 £j d + a d+a 

V £ / 2 > 2 



E ^ 2 " 2/t E 



m(z'), 



where n £ Z d_1 is a (d — l)-dimensional vector with 
integer components, and all the lengths are measured 
in units of lattice spacings. It is easy to show that 
max t exp(— ir 2 n 2 /t+t(z— z 1 ) 2 ) is attained at t = 7r|n|/|z— 
z'\ and is equal to exp(— 27r|n| \z — z'\). Because of this 
exponential decay in the last row of the above equation 
we will take into account only the terms with |n| = 1 and 
\z — z'\ = 1 (the corresponding improvements that take 
into account n = 2, 3, • • • and \z — z'\ = 2, 3, • • • are obvi- 
ous; as we will see even the contributions stemming from 
|n| = 1 and \z — z'\ = 1 are very small). It is clear that 
size dependent contributions that are due to the terms 
in the last row of i|A6|l will be exponentially small in L. 

For d = a = 3 the corresponding Mittag-Leffler func- 
tion can be expressed in the following simple form 



Taking into account that, if x > 0, 



W tE 3 . 3 (-t 3 )e- tx dt^g 3 (x), (A8) 
Jo 



where 



G 3 (x) 



v/3 arctan 



73 



In 1 



2x- 1 

1 1 
- + - 

X X 1 



In 1 



(A9) 



1 

it 2 



2e*/ 2 < 




we arrive at the following equation for the magnetization 
. (A7) profile 



arctanh \m*{z)\ = h(z) + K { c 2 m*(z) + c% n [m*(z + 1) + m*(z - 1)] + J2 G 3 (\z - z'\ 2 )m*(z') } , (A10) 



z f =0 
:' — jsl>2 



where 



C2 



T — 

1 + |n| ( 



3.602, 



(All) 



As we see, the contribution due to the first part in 
Eq. (|A1 2|> . and therefore the contributions due to the last 
row of Eq. 

It is easy to verify that 



6j are of the order of 1% in the constant cj™. 



and 



-— 7T 

3 



(-l) 1/3 X (V // 2-2zv / 37r) 

-(-l) 2 / 3 K (\f2 + 2iV3ir) 
0.00955 + 1.17354 « 1.183. 



•&(1) 



+ 3b- ln2 



(A12) 



G(x) ~ '-x- 2 - ^x- 5 + + 0(z- n ), x^oo. 

Zoo 

(A13) 

Setting Q = in Eq. (|A"To)) one obtains an equation 
having a form that is familiar in the mean-field theory 
of short-range systems. Actually in our case the system 
in question has short-range interactions in z direction 
and long-range ones within the planes perpendicular to 
z. The standard Ginzburg-Landau equation follows, for 
small m, after taking into account that arctanh(m) ~ 
m + m 3 /3 + 0(m 5 ) . A continuum version of the equation 
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FIG. 11: The behavior of the C?3-function as defined in <A9> . 



follows from the replacement m(z — 1) + m(z + 1) — > 
2m(z) + m"[z\. Obviously such a continuum version can 
also be constructed for the long-range system by adding 
the terms contributed by the function Q(x), which is, in 
this case, a continuous function. Note that the function 
Q (x) is well defined everywhere for x > and not only 
for x > 1 as we actually need it in the lattice formulation 
of the theory. Thus, in the continuum formulation of the 
theory the integration can be extended over the region 
z £ [0, L\. This does not change the long-range behavior 
of the magnetization profiles. In the continuum case the 
equation for the magnetization profile reads 



m*[z] + ^{m*[z)Y = h[z] + K i c 2 m*(z) + c%' 



2m* [z] 



dz 2 



'(z')dz' 



(A14) 



APPENDIX B: HOW TO SCALE 



In the discussion up to now we have tacitly assumed 
that the value of L is precisely known. This is, however, 
not only an experimental problem — due to the roughness 
of the surface, the existence of impurities, dust, etc. — but 
also a theoretical problem that might play a role when 
L is not "large enough". Let us make a brief comment 
about this issue. The definition of the size of the system 
is unambiguous only for systems with periodic boundary 
conditions. If iV is the number of layers with independent 
degrees of freedom, then the size of the system is simply 
L = Na, where a is the distance between the layers. 
Any point in the system is equivalent to any other point. 
Therefore, any layer is equally suitable to be taken as the 
origin with respect to which one measures distance. How- 
ever, how one proceeds for a system with (+, +) bound- 
ary conditions — when the first and the last layers have 
fixed degrees of freedom — is less clear. For consistency 
with periodic boundary conditions one can, of course, 
count the number of layers with independent degrees of 
freedom and let this, as in the case of periodic boundary 
conditions, be equal to N. Then the question is, shall we 
include in the total size of the system the half distance 
between the two outermost layers with independent de- 
grees of freedom and the adjacent layers with the fixed 
degrees of freedom? A reasonable approach seems to be 
that at least half of these distances should be taken to 
belong to the system, i.e. L = (N + l)a. This, of course, 
is not an unambiguous procedure. We have two layers 



with fixed degrees of freedom, which do not belong to the 
"substrate" surrounding the system that strongly prefers 
the ordered phase of the system in the case of (+, +) 
boundary conditions. Thus one has to somehow decide 
which portion of the bounding layers are to be counted 
within the system. Much more complicated is the case 
of systems with long-range interactions. Then any parti- 
cle (atom, or molecule) of the system interacts with any 
other one from the substrate. How then does one define 
the size of the system or the borderline between the sub- 
strate and the system? Let as denote by L the "true" size 
of the system (which we do not know) , and by L a reason- 
able approximation of that size (say, by taking, as above 
for the short-range case, the size to be L = (N + l)a). 
The last implies that L = L + S, where 5 -C L. Figure 
El demonstrates the difficulties when studying the scal- 
ing in systems with subleading long-range interactions 
of the van der Waals type. The role of the long-range 
surface potentials, which are irrelevant in the renormal- 
ization group sense but for moderate values of L, i.e. for 
a«L« L CI it (see Eq. (|1.11[) 1. contribute to the leading 
behavior of the finite-size susceptibility is clearly seen. 
One can say that, for such values of L, the quantity x s 
is a sort of "dangerous" irrelevant variable — in the sense 
that, despite being irrelevant, one cannot neglect it when 
L < -L cr it. We further note that the greatest deviation of 
the curves for different L from each other is around the 
maximum value of the scaling functions. The lack of the 
data collapse is due to the fact that x s ~ h w ^ s /^/L is not 
the same for all the curves (see Eq. I|l.ll|l ). Definitely 
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FIG. 12: The behavior of the finite-size susceptibility nor- 
malized per L A// " for different fixed values of h w ,s and for 
L = 500. Within the mean-field approximation one has that 
A /v = 3. The values of h W:S for which results are presented 
are h w>a = 0.1, 0.18, 0.3, 0.7, 1.1, 1.6, 2.0, 2.3, 3.0 and 4.0. They 
model the role of different substrates that surround a given 
fluid film with thickness L. 



similar spreading of the data for the finite-size suscep- 
tibility are to be observed if h w , s is kept fixed while L 
changes — say h WjS — 4 as in case of 3 He or 4 He confined 
by Au plates and— one considers L = 3000, 1000, 500, 100 
(see Fig. I13|) . However, when x s is kept fixed for the 



when L is "large but small enough" the effect of S, i.e. of 
the fact that we do not know exactly the system size L 
is clearly evident. In the case of \ one has x — £ 7 ^[1 + 
<S(7/^)L _1 + 0(S 2 L- 2 )]X x (x t ,x h \x s ,x b ,x UJ ), where x t = 
a t tLV», x h = a h hL^», x s = h w ^ d+2 - v)/2 ' a , x b = 
bL 2 ~ r '~' J , x^j — a u L~ u . Note now that the expansion 
of the above expressions in terms of L will yield all pos- 
sible non-universal (proportional to L^ 1 ) corrections to 
the leading finite-size behavior of the susceptibility with 
the greatest deviation of the curves for different 6 oc- 
curring near the maximum value of X x . For the sake 
of precision let us also note that similar corrections will 
be produced if one takes into account the change from 
L to L in the variables Xt,Xh,x s ,ib and x u . Thus, only 
the leading finite-size behavior can be determined unam- 
biguously. All corrections will depend on the definition 
of L, i.e. on S. There is, nevertheless, still something 
that one can do in order to check that the behavior of 
the susceptibility for L = 100 is simply due to the above 
explained unambiguity in the definition of L. Note, that 
if we normalize x to its value at a given point within the 
critical region - say to the value xo at the bulk critical 
point (T = T c , H = 0) - then, whatever the definition of 
L is, the leading behavior of the resulting quantity will 
not depend on this definition. Explicitly, we have 
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The result from the application of the above procedure is 
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FIG. 13: The behavior of the finite-size susceptibility normal- 
ized per L A ^ for different fixed values of L with the surface- 
field amplitude kept fixed at h W:S — 4. Systems with thickness 
L = 3000, 1000, 500 and L = 100 are considered. The curves 
represent the behavior of the susceptibility of 3 He and 4 He 
films with different thickness when the surrounding surfaces 
are made of gold. 

same values of L considered before, the data collapse im- 
proves greatly, as shown in Fig. 1141 One observes that 
violations of scaling are now clearly detectable only for 
the smallest system size L = 100. Let us now recall that 



FIG. 14: The finite-size susceptibility for L = 3000, L = 1000, 
L = 500, L = 250, L = 100. The amplitude of the surface 
magnetic field is rescaled in such a way that h w , s /v~L = const 
for all values of L and, as in the experimental realization, 
h w . s = 4 for L = 3000. One observes that practically all the 
curves for L > 250 coincide with each other, i.e. the scaling is 
indeed valid. The curve with L — 100 differs from the others. 
Thus L — 100 is too small, and in this case the importance of 
S is demonstrated. 

shown in Fig. We observe that all the curves, includ- 
ing L = 100, now allow for data collapse and that only a 
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small deviation of the curves from each other is observed 
for very large values of the scaling variable x t , when the 
onset of the nonuniversal corrections to scaling due to 
the role of the fluid-fluid interaction (i.e. proportional to 
Xb) is expected to set in. Thus, despite ignorance of the 
precise value of L, we are able to determine the leading 
finite-size behavior of the susceptibility. 



APPENDIX C: EVALUATION OF THE SCALING 
FIELD PARAMETERS FOR HELIUM 

We imagine a simple model for helium ( 3 He or 4 He) 
in which atoms interact via a pair potential w l (r,r') = 
— 4j'(r, r'). We assume that the fluid is bounded by a 
substrate whose particles interact with the helium par- 
ticles with a pair potential w l ' s (r, r') = — 4J l < s (r, r'). 
Within the lattice gas model for any given configuration 
C of particles {pf,p l j}, i € S, j e C with C and S de- 
noting the region occupied by the helium and substrate 
particles, respectively, the energy of the fluid is given by 

* = E + \ E <M (ci) 
= -4 E *M -> E J IM> 

where p l j G {0, 1} and pf e {0, 1} denote the occupa- 
tion numbers for the fluid and substrate particles, re- 
spectively. Since only the part {pj} belonging to the 
fluid becomes critical around T c and the fluctuations of 
the particles {pj} belonging to the substrate are unim- 
portant here, one can replace the latter ones by their 
mean-field values. If the fluid is in contact with a par- 
ticle reservoir with a given (excess) chemical potential fi 



and temperature T, the partition function for the liquid 
is 



Z = 



E ex p 
E cx p 



(C2) 



; 1 E AM 



ies,je£ 



+2 ]T JlM + nY.Pj 

ijec jec 



where p\ = (pf). Since the solid substrate is only weakly 
influenced by its surface, we assume (pf) = p s = const., 
in S so that 



z = ^2 ex p 

{pi} 



+2(3 J ljP l iP. 



(C3) 
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FIG. 15: The finite-size susceptibility for L = 3000, L = 1000, 
L — 500, L = 250, and L — 100 normalized to its value 
at T — T c . The amplitude of the surface magnetic field is 
rescaled in such a way that h w , s /\fL — const for all values of 

L- one finds for V. J,.'?: 
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where Zj > 1 characterizes the distance of the particle 
Pj from the boundary with the half space Si/ 2 occupied 
by the substrate. We consider the fluid particles to be 
in the region < z < L, where L is the width of the 
film confined between the two surfaces. Therefore, the 
partition function is 



E ex p 

{ P >} 




2 E J IM 



(C5) 



i.e., the system is equivalent to one with a spatially vary- 
ing chemical potential pj = p — Vj acting on a particle pj 
at a distance Zj + 1, < Zj < L, from the left boundary 
surface and at a distance [L + 1 — Zj) from the right one 
where Vj = V(zj) is given by the superposition 



V(z) = v s [(z + 1)-^ + (L + 1 - z)- 



with 



_ 47r (d-l)/2. 1 I 2 



PsJ 1 



(C6) 



(C7) 



In the current article we choose such boundary condi- 
tions that p(0) = p(L) — 1, where p{z) = (p l j). This 
is known as (+, +) boundary conditions. The pressure 
p in the fluid as a function of {pj} and T follows from 
ftp = \C\~ 1 InZ, where \C\ is the number of lattice sites 
in the region C. The critical properties of this model can 
be directly derived from the known critical behavior of 
the corresponding magnetic system that one obtains un- 
der the transformation m,; = 2j3j— 1, where m, 6 { — 1,1}. 
One arrives at 



0P = \{3\£\- 1 



X>+ E J L 



Pf, (C8) 



where / is the free energy of the magnetic system 



-ftf = In exp 

{m} 



E'"^'-' • -, E K > 



(C9) 



where 



K id = 0Jl d , and //, i i/,, • Y, ■ ■ ( C1 °) 



The mean-field critical properties of the model (|C9|) are 
well known. The critical exponents are /3 = v = 1/2, 
7 = 1 and for the uniform system the critical point is at 
{h = 0, K^ 1 = J^. Ki j}. At the critical point < TOj > c — 
m c = 0. Thus, for the fluid system we derive that the 
corresponding critical point is at {K = K Cl fi = p c }, 
where p c — —2J2iez d w hile at this point < pi > c — 
p c = 1/2. With the help of fi c and p c the corresponding 
expressions for hj can be rewritten in the form 







(P - Pc) 



(^+ 4 E4M 



8 27r( d - 1 )/ 2 
f (M - Pc) + ft—^TTT^T ( - 7 ' ""/'•- 



ar(^) 



x [(*,• + 1)" CT + (I + 1 - ^) _<7 ] • (Cll) 
From Eq. (|C11(1 (see also [3l]]) one identifies that 

^ -2^ d - 1 )/ 2 1 2 j 



0(p s J^ - p c J l ). (C12) 



The equation of the magnetization profile Ij2.2|l also 
directly follows from Eq. (|C9|) ■ Denoting m* — < mi > 
one obtains 



tanh 



E A ' / /,r + b - 



(C13) 



Taking into account that p c = 1/2 one can rewrite m* 
in the form m* — Ipi — 1 = {pi — p c )/Pc and, thus Eq. 
'131 takes the form 



tanh 



pj 



(C14) 



In what follows we take the 3 He or 4 He atoms to be 
constrained by an Au substrate. Then, according to 
Refs. HH0 d s ~ -270 meV A 3 /rg, where r is 
the distance between the helium atom and the Au sur- 
face. We will assume that ro is the same as the dis- 
tance between the 3 He or 4 He atoms (but being different 
for 3 He and 4 He cases, respectively). It is clear that 
ro provides the scale of the length of the unite cell of 
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the lattice on which we consider the fluid embedded. 
An estimation of r can be obtained from some gen- 
eral data for 3 He or 4 He. The critical density of 3 He 
is p c ~ 0.01375 mol/cm 3 ~ 0.04145 g/cm 3 35], while 
for 4 He it is p c ~ 0.017399 mol/cm 3 ~ 0.0690 g/cm 3 
[3H l37j , wherefrom one easily derives that at the criti- 
cal point one has 8.28 x 10 27 particles/m 3 for 3 He and 
1.38 x 10 28 particles/m 3 for 4 He. This leads to the conclu- 
sion that the space "allocated" for one particle is of the 
order of 120.77 A 3 for 3 He and of the order of 72.55 A 3 
for 4 He, i.e. the size of one 3 He atom at the critical point 
is of the order of 4.9 Awhile for 4 He it is 4.2 A. Thus, for 
v s one has v s ~ -2.3 meV ~ -3.68 x 10~ 22 J for 3 He 
and v s ~ -3.6 meV ~ -5.82 x 10~ 22 J for 4 He. For T 
near the critical temperature T c = 3.3 K of 3 He |3J| one 
has k B T c ~ 4.55 x 10~ 23 J and thus (3 ~ f3 c ~ 2.2 x 10 22 
J" 1 with f3 c v s ~ -8.1. For 4 He T c = 5.2 K jH and 
therefore k B T c ~ 7.17 x 10~ 23 J, (3 ~ (3 C ~ 1.4 x 10 22 
J -1 , and j3 c v s ~ —8.1. Taking into account that the 
atomic weight of Au is 196.97 u, whereas its density is 
19.3 g/cm 3 and having in mind that the atomic weight 
of 3 He is 3 u and that the atomic weight of 4 He is 4 u 
(where u — 1.66 05 x 10 -27 kg is the atomic mass unit), 
it is easy to verify that the number density of Au is 7.1 
times larger than the number density of the 3 He and 5.7 
times larger than that of 4 He at the critical point of the 
respective fluid. Since within the mean-field theory, the 
number density of both 3 He and 4 He at their respec- 
tive bulk critical points is p = p c = 1/2, we obtain that 
p s ~ 3.55 for 3 He films and p s ~ 2.85 for 4 He films. As an 
estimate of J l ' s one immediately derives from Eq. (|C7|) 
(for d = a s = 3) the result that that J l ' s ~ 4.95 x 10" 23 J 
for 3 He and that J l > s ~ 9.75 x 10" 23 J for 4 He. Next, ne- 
glecting the contribution due to J , i.e. the interaction 
between the atoms of 3 He and also between the atoms of 
4 He, one finds that h W:S ~ — \j3 c v s ~ 4.05 both for 3 He 
and 4 He films, i.e. the surface field is, indeed, relatively 
large and cannot be neglected. 

Next, we justify the approximation made for J . 
Within the mean-field approximation we have /3 C J — 

0. 160. Thus, from the experimentally known value of 
Pc we conclude that J 1 ~ 7.3 x 10~ 24 J for 3 He and 
J 1 ~ 1.14 x 10" 23 J for 4 He. Note that these estimates 
are very close to those based on the general expectation 
that k B T c ~ J 1 , which leads to J 1 ~ 10~ 23 J. There- 
fore, J 1 '" /J 1 ~ 6.6 for 3 He and J l > a /J l ~ 8.5 for 4 Hc, 

1. e. the interactions of the atoms of 3 He and 4 He with 
the Au substrate are much stronger than the interactions 
between themselves. If, nevertheless, one insists on tak- 



ing these interactions into account a simple calculation 
shows that h WjS changes from /i WjS = 4.05 to /i w-s = 3.96 
for 3 He and from h w , s = 4.05 to /i WjS = 3.97 for 4 He. 
Summarizing, one can conclude that the surface field has 
almost the same value for both the 3 He and 4 He films 
bounded by Au surfaces, and thus one can predict that 
the finite-size behavior of their finite-size susceptibilities 
for a given fixed L will be practically indistinguishable 
for both fluids. 

We finish this Appendix by briefly commenting on the 
correlation length amplitudes for the correlation lengths 



±i*i-* 



and 



e i (T) = ^(T^T c ± ,/ l = 0)~^ ± |i| 



Uh) = Hootch -> 0) ~ i , h \h\- v '^. 



One can show that in the case of a van der Waals fluid- 
fluid potential the amplitude £q~ of the second moment 
correlation length is |3lJ 



it 



Furthermore one has 



2d 2 



r l+r d 



E 



r l+r d + ° 



1/2 



(C15) 



£ -=£ +/V2, and £ , h = £ +/^3. (C16) 

The numerical evaluation of the sum l|C15(l for d = a = 3 
in the case of a simple cubic lattice then gives £o~ = 
0.635 a. Taking into account that, as we derived above, 
a = 4.9 A for 3 He and a = 4.2 A for 4 He we obtain that 
$ = 3.11 A for 3 He and that £+ = 2.67 A for 4 He. Of 
course the procedure used to calculate the above numbers 
constitutes a very strong approximation and one shall not 
expect to reproduce the best known values of these quan- 
tities. Nevertheless, the comparison with the known data 
reported in the literature $ = 2.71 A [H H3 for 3 Hc, 
and £q" = 2.0 A[3!| for 4 He, shows that the above results 
are not too bad. One straightforward way to improve 
the above approximation is to consider the fluid imbed- 
ded not on a simple cubic but on a body centered cubic 
lattice which is probably "closer" to the reality since then 
the atoms are more closely pact. For such a lattice we 
obtain £+ = 0.574 a, and thus £+ = 2.811 A for 3 Hc, 
and £q = 2.409 A for 4 He. These results are indeed es- 
sentially close to the ones obtained by using much more 
elaborate methods El HI HI. 
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